Portable Implementation of a Quantum Thermal Bath for Molecular Dynamics Simulations

Recently, Dammak and coworkers (Phys. Rev. Lett. 103:190601, 2009) proposed that the quantum statistics of vibrations in condensed systems at low temperature could be simulated by running molecular dynamics simulations in the presence of a colored noise with an appropriate power spectral density. In the present contribution, we show how this method can be implemented in a flexible manner and at a low computational cost by synthesizing the corresponding noise ‘on the fly’. The proposed algorithm is tested for a simple harmonic chain as well as for a more realistic model of aluminium crystal. The energy and Debye-Waller factor are shown to be in good agreement with those obtained from harmonic approximations based on the phonon spectrum of the systems. The limitations of the method associated with anharmonic effects are also briefly discussed. Some perspectives for disordered materials and heat transfer are considered.